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The techniques of Pauli-Villars regularization and discrete light-cone quantization 
are combined to analyze Yukawa theory in a single-fermion truncation. A special 
form of the Lanczos algorithm is constructed for diagonalization of the indefinite- 
metric light-cone Hamiltonian. 



1 Introduction 

As a step toward development of a method for nonperturbative solution of four- 
dimensional quantum field theories, we consider a single-fermion truncation of 
Yukawa theory in discrete light-cone quantization (DLCQ).u To regulate the 
theory we introduce Pauli-Villars (PV) bosons with indefinite metricla into the 
Fock basis-jThis extends earlier work on model theories by Brodsky, Hiller and 
McCartorHu to a more physical situation. 

The importance of Pauli-Villars regularization stems from its ability to 
regulate the continuum theory with counterterms generated automatically, ex- 
cept for a trivial mass counterterm. The numerical method, in this case DLCQ, 
can then be applied to a finite theory. The bare parameters are fixed via phys- 
ical constraints and become functions of the PV masses and of any numerical 
parameters. The original theory is recovered in the following sequence of lim- 
its: infinite numerical resolution, infinite (momentum) volume, and infinite PV 
masses. 

The DLCQ formulation is based on periodic boundary conditions for bosons 
and antiperiodic conditions for fermions in a light-cone box — L < x~ = (t — 
z) < L, — Lj_ < x, y < Lj_. The light-cone 3-momentum p = (p + ,p±_) is then 
on a discrete grid specified by integers n = (n, n x ,n v ): p + = (E + p z ) = ^n, 
p_L = {j^n Xl j^-riy). The limit L — > oo can be exchanged for a limit in terms 
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of the integer resolutionU K = ^P + for total longitudinal momentum P + . 
Also x = p + /P + becomes n/K, with n odd for fermions and even for bosons. 
The light-cone Hamiltonian H^c = P + P~ is independent of L. 

Because each n is positive, DLCQ automatically limits the number of 
particles to no more than ~ K/2. The integers n x and n y range between 
limits associated with some maximum integer N±_ fixed by the invariant-mass 
cutoff mf + p 2 ^ < XiA 2 for each particle i. We then have a finite matrix 
representation where integrals are replaced by discrete sums 

J d P + J d 2 P± f( P + ,p ± ) ~ ^ (jP) ^2f(nn/L,ft±%/L±) . (1) 

This trapezoidal approximation is improved through the inclusion of weighting 
factors El that take into account distances between grid point locations and the 
integration boundaries set by the cutoff A 2 . In the following sections we discuss 
how these techniques can be applied to Yukawa theory. 



2 Yukawa theory 

CQ Ham 
lion, istl 

^M 2 + 6M 2 + (n ± ir/L ± ) 2 



The DLCQ Hajmiltonian for Yukawa theory, when truncated to include only 
one fermion, 
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where M is the fermion mass, fi = fj,o the physical boson mass, ^ the ith PV 
boson mass, e\ = — (A, i) the polarization vector for helicity A, and 



dim, a jmJ 



(3) 
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Modes with zero longitudinal momentum have not been included. Fcrmion 
self-induced inertia terms are also not included, because they cancel between 
PV-boson terms. A fermion mass counterterm has been included to remove 
shifts proportional to In fii j fi. 

The number of PV flavors is three. Their couplings are given by £i<7, where 

6 = \f\Ck\ and 



1+5^ = 0, /i 2 +£ci/i? = 0, ^a/x 4 2 ln(^/M 2 ) = 0. 



(4) 



The sign of Ci determines the norm. This arrangement is knownEm to produce 
the cancellations needed to regulate the theory and restore chiral invariance in 
the M = limit. 

To Hlc we must add an effective interaction, modeled on the missing 
fermion Z graph, to accomplish cancellation of an infrared singularity in the 
instantaneous fermion interaction. The singularity occurs when the longitu- 
dinal momentum of the instantaneous fermion approaches zero. The effective 



interaction is constructed from the .pair creation and annihilation terms in the 
Yukawa light-cone energy operatorB 
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combined with the denominator for the intermediate state 
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To complete the cancellation of the singularity in the kinematic regime of 
positive longitudinal fermion momentum, the instantaneous interaction is kept 
only if the corresponding crossed boson graph is permitted by the numerical 
cutoffs. 

The single-fermion eigenstate of the Hamiltonian is written 
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with normalization 



$2 • <$> a = 16tt 3 P + 5(P' - P) . 



(8) 



The wave functions ^>S™ describe a contribution with no physical bosons and 
rii bosons of the ith PV flavor. 

Mass renormalization is carried out by fixing the mass M 2 of the dressed 
single-fermion state. This is imposed by rearranging the eigenvalue problem 
Hlc^o = M 2 & a into an eigenvalue problem for SM 2 : 
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where the new wave functions are related to the originals by <j> = (f>/y/x and K. 
is the interaction kernel. 

To fix the coupling we set a value for the expectation value (:</> 2 (0):) = 
:0 2 (0): t&o- for the boson field operator (j>. From a numerical solution it can 
be computed fairly efficiently in a sum similar to a normalization sum 
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The constraint on ( 
eigenvalue problem. 



6 2 (0):) is satisfied by solving it simultaneously with the 



3 Diagonalization method 

The method used for diagonalization of the eigenvalue problem is based on 
a special form of the Lanczos algorithm U designed to handle the indefinite 
norm. Let r\ represent the metric signature, so that numerical dot products 
are written (</>'\<j>) =(/)'*■ rj<f>. An operator H is self-adjoint with respect to this 
metric if O H = = H . The Lanczos algorithm for the diagonalization 

of H then takes the form 

Qj+l = rj/Pj , fj = Hq - Jj-icfj-i - a 3 q 3 

Vj+i = sign(f/ ■ rfj) , v x = sign(g* ■ rjqi) , (11) 
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The original (large) matrix H acquires a new (small) matrix representation T 
with respect to the basis formed by the vectors qj : 



H ->T 
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The elements of T are real, and the matrix is self-adjoint with respect to 
the induced metric v. We can solve Tci = XiCi for eigenvalues and right 
eigenvectors and find that H(j)i ~ \i<pi, with 
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Extreme eigenvalues are well approximated after only a few iterations. 

Because of the indefinite metric, the physical one-fermion state is not nec- 
essarily the lowest mass state. It is instead identified by the following char- 
acteristics: a positive norm, a real eigenvalue, and the largest bare fermion 
probability between and 1. Each of these characteristics can be computed 
without constructing the full eigenvector, provided one saves the first compo- 
nent of each Lanczos vector qj to reconstruct the bare fermion probability. 



4 Future work 

Current work on this formulation of Yukawa theory is focused on the tuning 
of DLCQ weighting factors, both in general and specifically with respect to 
the infrared cancellations between instantaneous fermion interactions, crossed 
boson graphs, and the effective Z interaction. The singular interactions make 
the numerical representation rtupre sensitive to the weighting than was the case 
for earlier model calculations. 0a 

Within the present no-pair approximation, we can next consider two- 
fermion states. This would allow consideration of a true bound state.113 For the 
full theory, with fermion pair creation, we can again consider the one-fermion 
and two-fermion sectors. This will require further analysis of the infinities of 
the theory and determination of the appropriate PV particle types and inter- 
actions. 
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Other theories that could be considered include quantum electrodynamics, 
where one. might use PV regularization to repeat a calculation by Hiller and 
BrodskyEU of the electron's anomalous moment for large coupling. Quantum 
chromodynamics will require a somewhat different approach; a broken super- 
symmetric formulation, with massive partners playing the role of PV particles 
may be the correct route. 
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